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As shown by Gertsenshtein in 1961, an external magnetic field can catalyze the mixing of graviton and 
photon states in a manner analogous to neutrino-flavor oscillations. We first present a straightforward 
derivation of the mechanism by a method based on unpublished notes of Freeman Dyson. We next 
extend his method to include boundary conditions and retrieve the results of Boccaletti et al. from 


1970. We point out that, although the coupling between the graviton and photon states is extremely 
weak, the large magnetic fields around neutron stars ~10!4 G make the Gertsenshtein effect a plausible 
source of gravitons. We also point out that axion-photon mixing, a subject of active current research, 
is essentially the same process as the Gertsenshtein effect, and so the general mechanism may be of 
broad astrophysical and cosmological interest. 
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1. Introduction 


In a classic 1961 paper [1], Mikhail Gertsenshtein demon- 
strated that electromagnetic waves passing through an external 
magnetic field in curved space could be transmuted into gravi- 
tational waves. The result was entirely classical but in quantum 
mechanical language one can say that the external field “cat- 
alyzes” a resonant mixing between photon and graviton states 
in a manner analogous to the mixing of neutrino flavors. Al- 
though the coupling between the graviton and photon states is 
extremely weak, the large magnetic fields associated with, for 
example, neutron stars, make the Gertsenshtein effect a plausible 
source for producing gravitons [2]. Zel'dovich [3], however, and 
Dyson [4,5] pointed out that the Gertsenshtein mechanism re- 
quires coherence between the gravitational and electromagnetic 
waves. Large magnetic fields create electron—-positron pairs from 
the vacuum, which in turn alters its index of refraction. The speed 
of electromagnetic wave propagation is lowered relative to that 
of the gravitational wave and the Gertsenshtein mechanism is 
quenched. 

The Gertsenshtein effect received some attention in the 1970s 
when Boccaletti et al. [6] provided a comprehensive discussion 
and de Logi and Mikelson [7] calculated the cross section for 
graviton-photon conversion. Cross sections have also been com- 
puted by Long et al. [8], and the effect has been discussed by 
Ejlli [9] as well as Jones and Singleton [10]. More recently, the 
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mechanism has been mentioned in the context of the primor- 
dial gravitational wave, or graviton, background by Vagnozzi and 
Loeb [11] and various others [12-16]. Nevertheless, most physi- 
cists apparently remain unaware of the Gertsenshtein effect, al- 
though current calculations of axion-photon conversion employ 
a nearly identical formalism (see Section 4). 

As it turns out, Freeman Dyson rediscovered the Gertsenshtein 
effect in 2005 [4], evidently unaware of the previous work. When 
informed that he was not first, he chose not to publish his notes, 
but gave them to one of us (TR). As Dyson’s derivation is simpler 
than the others we have seen and has remained unpublished, 
we thought it worthwhile to make it available in some form. In 
the following section of this paper we present Dyson’s approach, 
filling in some steps and clarifying a few ambiguities. In Section 3 
we make contact with Boccaletti et al.’s results. It will be seen that 
Dyson’s method is much simpler. Finally in Section 4 we point 
out the similarity to axion-photon conversion, and conclude by 
indicating some possible applications of the Gertsenshtein effect 
in the cosmological context. 


2. Graviton-photon oscillations 


Consider the propagation of electromagnetic and gravitational 
waves in a flat spacetime. For the gravitational wave (GW), we 
assume as usual that it is represented by a weak perturbation 
hyv(x, t) < 1 traveling on a Minkowski background such that 
the full spacetime metric is 


Zu = Nuw thw ; g =a" h ; h” = nè” hap, (2.1) 


where we take the flat-space metric to be n,,, = (—1, 1, 1, 1). 
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We carry out our calculations in the standard transverse- 
traceless (TT), or Lorentz, gauge. If the wave is taken to be 
propagating in the z-direction, then h,, = h,,(z,t) and the 
transverse components of the field are hy, = —hg2 and hy2 = hy, 
which represent the two independent polarizations of the GW.! 

In the Lorentz gauge, the linearized gravitational wave equa- 
tion (see [17], chap.18 and [18], chap.10) is 


hyy = —167 Tv, (2.2) 


with O = “ð ,. We assume that the energy-momentum tensor 
Ta Consists entirely of the electromagnetic energy of the system. 
In doing so, we ignore gravitational backreaction, which is a good 
approximation as long as h,,, < 1. Therefore we have 


1 
AnT yy = FuoFy* — ggo FaF”, (2.3) 


where F,, = iAy = a is the usual electromagnetic field tensor 
and A,, = (—¢, A) is the vector potential. In terms of the electric 
and magnetic fields, the Maxwell stress tensor reads 


1 
4r Tij = —(E;iEj + BiB;) 4 5 SHE -E + B-B). (2.4) 
With our conventions, Fj = Aji — Aij = £ijkBę and Fio = Ao,i — 
Aio = —d¢ — Ai = E,.2 (Here, Ejx is the usual permutation 
operator.) 


We further assume that the electromagnetic wave obeys the 
vacuum Einstein-Maxwell equations: 

1 
——( VEF"), p = 0. 
VE i 
Since with metric (2.1) the determinant is g = 1 + O(h;) x 1, 
(2.5) amounts to 


(2.5) 


(gg Fus) g = 0. (2.6) 
Working to first order in h, we have 
n“ n” Favs g Me ne Fay po Fuh”, p 

— PhO Fav, g PF yyh" 9 = 0. (2.7) 


However, in the Lorentz gauge hô”, p= 0, so the third term 
vanishes. Discarding the second and fourth terms as small, but 
keeping the final term as potentially large leaves 


age Fans —nP’ hg Fuy = 0. (2.8) 


For a = 0 (2.8) gives merely V - E = 0, as in flat space. For a = i, 
the equation becomes 


—E; + eixeBesx +hiyEj — hijsk &jneBe = 0. (2.9) 


(Since hË = hj, it is immaterial whether one writes this equation 
with indices up or down.) The È; term can be eliminated as in 
freshman physics to get a wave equation for B. Taking the curl of 
(2.9) yields 


—8piEisp +E rpi€ike Bes kp +e rpi(hiyE;),p —ErpiEjke (hij, k By )sp = 0. (2.10) 
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B; tt. By virtue of the homogeneous Maxwell equation V - B = 0, 
the second term is V x V x B= —V°B = —B, gv 

Next, we assume that the electric field E is entirely that of the 
electromagnetic wave, but that the magnetic field consists of the 
field b of the electromagnetic wave plus a background field Bo, 
such that B = Bo + b with Bo > b; Bo > E. Thus we can ignore 


the third term above as second-order small and (2.10) becomes 
Br, tr —Br sce —Expi€je(hy sk Be),p = 0. (2.11) 


With the same assumptions as above and discarding terms 
~ b?, the gravitational wave equation (2.2) becomes 


Oh; = 4[Bo:Bo; + Boibj + biBoj — Tae, + 2Boxbx)I. (2.12) 

The trick is, first, to choose Bọ = constant,* which linearizes 
both (2.11) and (2.12), and secondly, to choose Bo at right angles 
to b, which ensures a nonzero coupling for the h42 mode in (2.12). 
Assuming that the electromagnetic wave is traveling in the z- 
direction, we take b = b,(z, t) and Bo = Box. Then from (2.11): 


by — by = h'zBo, (2.13) 


where the prime (’) indicates differentiation with respect to z, and 
from (2.12): 


hy — hi, = —A4Boby. (2.14) 


Note that the background field Bọ makes the mixing between 
a vector field (photon) and a tensor field (graviton) possible by 
providing the missing angular momentum. Given that Ohj ~ 
Ty ~ E,E; + B;B;, for the given field configuration the only possible 
linear mixing is between the h,2 component of the gravitational 
wave and the b-field of the electromagnetic wave. There is no 
interaction with the other polarization state hyo 

It remains now only to solve the two coupled equations (2.13) 
and (2.14). We search for “wave packet” solutions f the form 
f(z, t) =f-(z —t)f,(z + t), where r stands for “rapid” and s stands 
for “slow”. The Ur-example would be 


f = 
= 5 [cos (ky + k;)z + (@s — wr )t) + cos ((k; — ks)z 
— (a, + @s)t)], 


cos (k;z — œt) cos (k,Z + wst) 


where œ >> œ; and in our units k, = œr; ks = w;. Such wave 
packets can thus be regarded either as product functions con- 
sisting of a high and low frequency component or as a coherent 
superposition of waves traveling in the same direction, which can 
be expected to form beats as in elementary physics. 

Dropping subscripts, we write h = h,(z — t)h,(z + t) and 
b = b,(z — t)bs(z + t), where in our units h. = —h, and h, = hs, 
with similar expressions for b, and b,. Inserting the trial functions 
into (2.13) and (2.14) gives 


B B 
bib, = z (hh + QW + hh!) ~ -7h (2.15) 


By virtue of the homogeneous Maxwell equation V x E = B, 
which is the same as in flat space, the first term is 3?B,/ðt? = 


1 We use units in which G = c = 1 and follow MTW [17] conventions 
throughout. In particular, Greek indices = 0...3 are spacetime indices, while 
Latin indices are spatial indices = 1...3. Repeated indices in any position are 
summed. 

2 Because hy, = —h22, we must also have T;; = —T22, which is not generally 
true for the stress tensor (2.4) unless E; = B3 = 0. This will be the case 
for transverse electromagnetic waves traveling in the z-direction (“radiation 
gauge”), which we also assume. 


3 Here we have used the vector identity V x V x A= V(V-A)— V2A. 

4 Note that for i = j = 3, (2.12) gives Oh33 ~ Bi # 0. This entails the 
appearance of a longitudinal component of the gravitational wave h33, induced 
by T33. Given that in our case Oh33 ~ Bi = constant, this mode does not mix 
with the other two transverse components; therefore, for the purpose of our 
derivation, we can effectively ignore it. 

5 Another way to see this is to note that the polarization tensors for the GW 
can be written as the dyadics e, = ii —ĵj and e, = ij + ji. Thus, e, - Bo || Bo 
and e,-b || b, so the hy, mode produces no interaction between By and b, 
whereas ex : By || b and e, -b || Bo, so the hj. mode does allow interaction 
between Bo and b. 
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and 


hh, = Bob; bs. (2.16) 


In (2.15) we have approximated h” as hi/hs, since by assumption 
h, is rapidly varying and h, is slowly varying. 
The equations can now be solved by inspection. Pick 


3/2 
B B 
kK. = An ; k, = ay/Bob, ; b; = -f hs, (2.17) 


with a an arbitrary constant. These choices satisfy (2.15) and 

(2.16). The functions b, and h, which describe the shape of the 

wave-packets can also be chosen arbitrarily and have no effect on 

the behavior of the functions b, and h,, which describe the slow 

oscillation of the packets between photon and graviton states. 
Letting a = 1/2 gives 


, 1 1 
hl = 5 V Bobs ; b= = 5B hs, (2.18) 
which have solutions 
t t 
h; = Asin =) ; bs = Ay Bo cos =) : 
(2.19) 
provided that L = 2/Bpo. Notice also that 
b2 
h? + A = A’ = constant, (2.20) 
0 


which shows that the total energy remains constant as the packet 
oscillates (see the Appendix for details). The frequency of oscilla- 
tion is given by the magnitude of the magnetic field: œs = L7! = 
Bo/2. 

We have seen that a uniform classical magnetic field Bo cat- 
alyzes a linear mixing of the gravitational and photon fields. A 
quantum of the mixed field oscillates between photon and gravi- 
ton states with a mixing length L independent of wavelength, 
namely 


2 

= 3 
where Bo is the component of the background field perpendicular 
to the direction of propagation of the quantum. This means that 


if a single photon travels a distance D through the uniform field 
Bo, it will emerge as a graviton with probability 


P =sin’(D/L). 


L (2.21) 


(2.22) 


The quadratic dependence of P on D makes this process interest- 
ing as a possible astrophysical source of gravitons. 
In cgs units 


2c? 1 Gauss 
L= x 2Mpc : 
Bo 


(2.23) 


Around a magnetar, the field can be ~10"4 G, leading to a mixing 
length of only ~10° km. From (2.22) above, an observer stationed 
within the boundaries of such a field at a distance D >> L from 
the star could in principle observe a significant periodic change of 
magnetar luminosity on a timescale of seconds. More realistically, 
the extent of a neutron star’s magnetic field is ~10 km. Beyond its 
field boundaries oscillations cease (Section 3) and the brightness 
of the star is diminished by a factor of 10~!° in the direction 
perpendicular to the magnetic field. Conceivably, two observers 
could infer the presence of the Gertsenshtein mechanism by 
comparing photon fluxes in directions parallel and perpendicular 
to the neutron star B-field. However, since the effect is classical, 
such an observation would not constitute a direct detection of 
gravitons [2,19,20]. 
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We also note that 10 G is approximately the Schwinger 
limit, at which one expects the electromagnetic field equations 
to become nonlinear and electron-positron pairs to be created. 
Pair creation would in turn alter the index of refraction of the 
vacuum and change the speed of light, destroying the presumed 
coherence of the gravitational and electromagnetic waves. Hence, 
for such strong fields one should do a more detailed calculation 
before drawing definite conclusions. 


3. Inclusion of boundary conditions 


In the previous section we obtained a general solution for 
graviton-photon oscillations. Boccaletti et al. [6] employed a 
more conventional equation-solving method to obtain a solution 
in which conversion takes place within a finite region, say for 
0 < z < D. Such a scenario is probably more relevant to graviton 
production in the vicinity of neutron stars, for example. 

We now employ the Dyson approach to retrieve their results 
for photon-to-graviton conversion in this situation. 

We assume that the constant background magnetic field Bp = 
Box is confined to the finite region 0 < z < D (region II), with 
a plane right-moving electromagnetic wave traveling along the 
z-axis incident upon it. By the Gertsenshtein effect, the electro- 
magnetic wave produces a left-traveling gravitational wave in 
region! (z < 0), and a right-traveling gravitational wave in region 
Ill (z > D). 

In the region where the magnetic field is non-zero, the trans- 
verse component of the gravitational wave h42 is described by a 
wave packet of the form,° 


hj(z + t)h,(z — t) = Asin (@,(z + t)) sin (@,(z — t) (3.1) 


where œr (@; = L7!) is the frequency of the rapidly (slowly) oscil- 
lating component of the metric, A is the maximum amplitude of 
the wave, and for simplicity we have set the global phase ¢ = 0 
in (2.19). 

The most general solution for the gravitational wave in region 
II is then a superposition of wave packets of the form (3.1) and 
plane waves. Following [6], we turn now to a complex formalism 
to describe these waves, with the understanding that physical 
solutions are always given by the real part of complex amplitudes. 

Using this formalism, h,2 in region II is given by (the real part 
of) 


WL, t) = qe@r(z-t) + be ert) + Ae i@s(z+t) el@r(z—t) 

for 0<z<D. (3.2) 
The first two terms represent right-traveling and left-traveling 
plane waves respectively, while the third term is a superposition 


of wave packets.’ 
In both regions z < O and z > D where the background 


magnetic field is zero, the wave equation for the metric is Oh, = 
0, and the corresponding vacuum solutions are 

hl (z, t) = ce™®rC+O for 20, (3.3) 
and 

h! (z, t) = de®- for z >D, (3.4) 


describing a generic left-traveling (z < 0) or right-traveling (z > 
D) plane wave. 


6 As we saw in Section 2 the function h,(z — t) can be chosen arbitrarily. 
Here we make the simplifying assumption that h,(z_;) « sin(@,(z — t)). 

7 Specifically, the third term is a combination of the four elemen- 
tary wave packets cos(w,(z + t))cos(w,(z — t)), cos(w,(z +t))sin(@,(z — t)), 
sin (ws(Zz + t)) cos (w (z — t)), and sin(w,(z + t)) sin(@,(z — t)). 
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In order to determine the constant coefficients (a, b,c, d) in 
terms of A we impose the usual boundary conditions demanding 
continuity of the functions and their derivatives at z = O and 
z=D: 


h!(0, t) = hi, (0, t), 
azh; (0, t) = ðzh (0, t), (35) 
h(D, t) = hi} (D, t), 
a,hi,(D, t) = a,hi}(D, t). 
Or, in terms of the coefficients: 
c=a+b+ Ae, 
c=b-a (: =) Ae lest 
Sr (3.6) 


d=a+ be72i@rD + Ae ies(D+t) 
d= a—ben2irD p qeit (4 _ Os 
Or 


This is a system of four linear equations for five unknowns. In 
terms of A, coefficients (a, b, c, d) are: 


Ws L; 
a=—A(1— east 
( oa) 


ho AS. 


eiD(2or — 0s) p—iwst : 


20, 


Os iD(2wr—ws)\ p—iwst 
c=A 1 — e ores] eT lost 
Jo ) 


d= a(t Z Ja e7 iPos )p-iost, 
r 


(3.7) 


2w 


The presence of the time-dependent factor exp(—iw,; ) means that 
the coefficients in (3.7) are not truly constant. However, given our 
assumption that e~s‘ is slowly varying, to first approximation 
we can consider them time-independent: 


a=—Aj|1 Ps 
~ 2a, )’ 


b= — AOS ePRoros) 
20, 


3.8 
Os a > eiD(2er—ws)). ( ) 


d= a(i æ )q en iDos 
20, 


One can easily check that the metric hj with the coefficients in 
(3.8) satisfies the wave equation in all three regions (in our ap- 
proximation) as well as the boundary conditions in (3.5). Hence, 
the full solution for h42(z, t) is 


(62) , j 
aV = elD(20r—Ws))p—tor(Z+t) 
r 


hi(z,ty=A 


h! (z, t) =—Aļ1— as ert) n A— ePlor-os)o7ior(z+t) 
20, 20; 


AeH) pior (z—t) 


, 


Os Ja e7 Dos pior(2—t)_ 
' 


20 


hi(z,t)=—A (1 
(3.9) 


We can similarly derive a solution for the electromagnetic 
wave. From Eq. (2.17), by = hi. /(2./Bo) and b; = 2h! /~/Bo. Given 
that h,(z—t) = el@r@- and h,(z+t) = Ae!@s@+9, we have b,(z— 
t) = iw, /(2/Bo)e" 2-9 and b.(z + t) = —2iAws//Boe 2+, so 
the electromagnetic wave packet in region II is of the form 


1 ; . 
b,(z — t)b,(z + t) = sre merle, (3.10) 
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By assumption, the incident electromagnetic wave is right- 
moving in both regions I and III. In addition, we need to include 
a produced left-moving electromagnetic wave in region I, as well 
as a combination of right-moving and left-moving plane waves 
in region II, as we did for the gravitational wave in (3.2). Con- 
sequently, after imposing the relevant boundary conditions as in 
(3.5) and disregarding the slowly oscillating term exp(—iw,t) % 1, 
the full solution for b(z, t) is: 


1 


I = _pior(z—t) L (1 — piD(2or—ws) —iwr (z+t) 
b (z, t) = Bære + z” e )Aase : 


b"(z, t) E (i a )] werz- _ peter Aangerioret 
Or 


4 Aopen tasara), 


b"(z, t) E s(t ws Ja z wp, eier(—t), 
T 


2%, 


(3.11) 


where B is the incident electromagnetic wave amplitude. Thus 
far, the coefficient B is undetermined because we have had to 
solve for five coefficients but have only four boundary conditions. 
However, one can easily check that the wave equation (2.14) for 
hy2 in region II is consistent with the second of (3.11) only if 
B = A/2. One can also check that the resulting expressions for 
hı2 and b conserve energy. 

The solutions (3.9) and (3.11) describe the conversion of an 
electromagnetic wave to a gravitational wave, catalyzed by the 
constant background magnetic field Bo. In region II where By Æ 0, 
the electromagnetic and gravitational amplitudes slowly oscillate 
between each other via the conversion term exp(—iaw,(z + t)). 
Outside, in regions I and II, the oscillation stops and the wave 
packets propagate independently. 

The solution of [6] is an approximation of ours for ws < wr. 
In fact, to leading order in w,, the gravitational wave solution in 
(3.9) becomes 


w j A 
hl (z, t) = A5 S 1 _ e2iDor je tor(Z+t) 
Or 


hy(Z, t) = A w, elor(z—t) A e2iDor e—iwr(z+t) iw ZA ior (Z—t) 
2 S 
Or 20; 


hZ, t) = —iDo, Ae), 
(3.12) 


where we have used e~s? ~ 1 — iw,z for small as. 

The coefficients of (3.12) are the same as in [6], given that 
the authors use a convention in which the wave packet solution 
is given by Az/(2iw, )exp(iw,(z — t)), with à = 2w,w,A. In this 
approximation, the amplitude of the gravitational wave in region 
I vanishes if the condition exp(2iDw,) = 1 is satisfied, and the 
amplitude of the gravitational wave in region III depends linearly 
upon the length and magnitude of the static field. 

Similarly, the approximate solution for the electromagnetic 
wave is, to leading order in ws, 


1 P 1 ; : 
b! (z, t) = 5 tween + ai! = e2iDar J Awe ret) 
1 : i or ; 1 
b! (z, t) = -Aw eE- — —¢2iPor Ay e EH) — i Ao wget’, 


1 ; 
b"'(z, t)= 5 (1 — iDows) Aw, eE», 
(3.13) 


If the condition e?”P®r = 1 is satisfied, the back-propagating wave 
amplitude vanishes and only the incident wave is left in region 
I. The solution in this approximation describes an incident elec- 
tromagnetic wave of amplitude A and energy œr, that propagates 


A. Palessandro and T. Rothman 


into region II and slowly converts to a gravitational wave with 
a rate given by œ = L~!. When the wave enters region III, the 
incident electromagnetic wave amplitude is effectively reduced 
by a factor of iDw, while the gravitational wave amplitude is 
increased by that same factor: 


D 
EM : A> a(1-i2), 


" (3.14) 
GW: 0> -iA 


4. Axion-photon conversion 


Quite generally the Lagrangian density for the interaction of 
gravitational waves with matter is 
L= -Sht (4.1) 
where in our problem T”” is given by (2.3) and (2.4). 

Currently, there is a large interest in axion—photon conver- 
sion [21-23], as axions are presently the leading candidates for 
dark matter. In a process much like the Gertsenshtein mecha- 
nism, axions emerging from the sun can be transformed into 
X-ray photons by a strong laboratory magnetic field and subse- 
quently detected by an X-ray telescope. Indeed, the axion-photon 
interaction Lagrangian can be written as 


Say 


E= <P Fw P, (4.2) 


where F“” = 1 uvoB Fag is the dual electromagnetic field tensor, 
a is the axion amplitude and ga, is the axion-photon coupling 
constant. 

In terms of the electric and magnetic fields, the above ex- 
pression amounts to £ = ga, AE B, which is essentially the 
same Lagrangian as the one we have been using except that the 
external magnetic field couples to the E-field of the photon rather 
than the B-field. 

Not surprisingly, then, we can employ Dyson’s approach for 
axion-photon mixing. In the high-energy limit in which one can 
ignore the mass of the axion, and in the presence of a background 
magnetic field Bp, the coupled equations for the axion a and the 
E-field of the photon e are analogous to (2.13) and (2.14): 


N 


ü — a" = —4ga, Boe 


4.3 
è = e" = Bay" Bo, ( ) 
and they admit oscillatory solutions with mixing length 
2 
L= ; (4.4) 
Zay Bo 


which agrees with the high-energy limit of [21]. As one can see, 
this is for all intents and purposes the result of Section 2. 


5. Conclusions 


The Gertsenshtein effect describes the conversion of electro- 
magnetic waves into gravitational waves and vice-versa. Conver- 
sion is made possible by an external magnetic field that provides 
the extra angular momentum necessary for a spin-1 field to mix 
with a spin-2 field. 

In this paper, we have presented a novel approach to deriving 
the Gertsenshtein effect due to Freeman Dyson. The key step 
in his method is to decompose the simultaneous solution of 
the linearized gravitational wave and Einstein-Maxwell equa- 
tions into rapidly and slowly oscillating components. The rapidly 
oscillating component describes the “shape” of the electromag- 
netic and gravitational wave packets, while the slowly oscillating 
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component describes the “beats” between the two fields, with 
the resulting oscillations having a mixing length independent of 
wavelength. 

In Section 3 we extended Dyson’s method to include boundary 
conditions and retrieved the result of [6] as a special case of 
our solution. We then employed (Section 4) the same method in 
the context of axion-photon conversion with essentially identical 
results. Given the quadratic form of Lagrangians in general and 
the similarity of the axion-photon Lagrangian to the gravita- 
tional interaction Lagrangian in particular, we may conjecture 
that Gertsenshtein-like mechanisms which cause mixing among 
multiple fields are “universal” and have widespread cosmological 
utility. Our approach can almost certainly be applied to any 
scenario involving these mechanisms. 

In future papers we plan to extend the effect to the (orig- 
inal) Yang-Mills SU(2) field as well as consider the possibility 
that graviton-photon mixing in the long-wavelength limit might 
generate something resembling a cosmological constant. 
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Appendix. Conservation of energy 


In this appendix we demonstrate in detail the statement made 
in Section 2, that energy is conserved during the oscillations 
between the electromagnetic and gravitational fields. 

Choosing for simplicity h,(z—t) = sin(w,(z — t)) as the rapidly 
oscillating component of the metric, and retaining factors of G, 
the gravitational wave packet is given by 


h(z, t) = h.(z+t)h,(z—t) = Asin (œw;(z + t))sin(@,(z — t)), (A.1) 


while, by virtue of (2.17), the electromagnetic wave packet is 
given by 


b(z, t) = b,(z + t)b,(z — t) = Awr cos (@s(Z + t)) 


1 
2G 
x cos(œ@r(z — t)). (A.2) 


Thus the slowly oscillating components are out of phase by x /2, 
but the rapidly oscillating components are as well, and so the 
above equations do not in fact produce beats. The reason is that 
(A.1) is not the physical gravitational wave packet. 

To see this, note that the energy density of the electromagnetic 
wave is 


b? 1 22 one? 2 
EEM = e 32G“ w; COS*(ws(Z + t)) cos*(w,(z — t)), (A.3) 
while the energy density of the gravitational wave is 
h2 


x EP sin?(w,(z + t))cos?(w,(z — t)), (A.4) 
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where we have approximated h =% Aw sin(@,(z + t)) cos 
(@r(z — t)) given that by assumption w, « œr. Eqs. (A.3) and 
(A.4) do describe coherent oscillations in the energy density of 
the two fields, suggesting that we take for the physical solution 
of the gravitational wave packet 


h 1 
2/G 2VG 
Indeed, the wave packets (A.2) and (A.5) describe a right-traveling 
plane wave of frequency œr that slowly oscillates between gravi- 
tational and electromagnetic states with frequency ws. Summing 
(A.3) and (A.4) and averaging over time gives 
_ Aw? 

~ 647G’ 


showing that the total energy of the system is conserved as the 
packet oscillates. 


h= Aw, sin(w,;(z + t)) cos (w,(z — t)). (A.5) 


(A.6) 


1 2 
(eem + ecw) = — (b° +h’) 
8x 
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